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Abstract 

Tadpole improved Wilson quark actions with clover terms on anisotropic lattices 
are studied numerically. Using asymmetric lattice volumes, the pseudo-scalar meson 
dispersion relations are measured for 8 lowest lattice momentum modes with quark 
mass values ranging from the strange to the charm quark with various values of the 
gauge coupling (3 and 3 different values of the bare speed of light parameter v. These 
results can be utilized to extrapolate or interpolate to obtain the optimal value for 
the bare speed of light parameter v ovt {ra) at a given gauge coupling for all bare 
quark mass values m. In particular, the optimal values of v at the physical strange 
and charm quark mass are given for various gauge couplings. The lattice action 
with these optimized parameters can then be used to study physical properties of 
hadrons involving either light or heavy quarks. 

Key words: Non-perturbative renormalization, improved actions, anisotropic 
lattice. 

PACS: 12.38.Gc, 11.15.Ha 



1 Introduction 



It has become clear that anisotropic lattices and improved lattice actions are 
the ideal candidates for lattice QCD calculations involving heavy objects like 
the glueballs, one meson states with non-zero three momenta and multi-meson 
states with or without three momenta. It is also a good workplace for the 
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study of hadrons with heavy quarks. In this work we present our numerical 
study on the quark action parameters suitable for heavy flavor physics. The 
gauge action employed in this paper is the tadpole improved gluonic action 
on asymmetric lattices: 

BTrPij 1 TrRij 1 TrR^' 
9l^f " 36 iul ~ 36 iul 

9 u 2 s 36 u A s \ U 

where P^- is the usual plaqette variable and R^j is the 2x1 Wilson loop on the 
lattice. The parameter u s , which we take to be the forth root of the average 
spatial plaquette value, incorporates the so-called tadpole improvement and 
£ designates the (bare) aspect ratio of the anisotropic lattice. With the tad- 
pole improvement in place, the bare anisotropy parameter £ suffers only small 
renormalization which we neglect in this study. Using this action, glueball and 
light hadron spectrum has been studied within the quenched approximation 
[1,2,3,4,5,6,7]. 

It has been suggested that relativistic heavy quarks can also be treated with 
the help of anisotropic lattices (the Fermi lab approach), possibly with im- 
provements [8,9,10,11,12]. Using various versions of the quark actions, charmed 
meson spectrum, charmonium spectrum, charmed baryons have been studied 
on the lattice [13,14,15,16,17,18]. Another type of application of the anisotropic 
lattices is the calculation of hadron-hadron scattering lengths within the quenched 
approximation [19,20,21,22]. However, in order to take full advantage of the 
improved quark action on anisotropic lattices, some parameters in the action 
have to be determined, either perturbatively or non-perturbatively, in order to 
gain as much improvement as possible. Some numerical studies of these param- 
eters have already appeared in the literature [23,12,24,25,26]. The anisotropic 
quark actions used in these studies fall into two categories. These two cases 
differ mainly in the choice of spatial Wilson parameter r s . According to the 
tree-level study [8,11], the choice of r s = l/£ has a virtue that the optimal 
parameters in the action as a function of the quark mass contains no correc- 
tions of the form moa s . The quark mass dependent corrections comes in only 
in terms of moat which is assumed to be small. As a result, the optimal values 
of the parameters can be approximated by their values in the zero quark mass 
limit. That is to say, tuning of the parameters in the action becomes almost 
quark mass independent. The disadvantage of this choice is that the doubler 
are not very well separated from the ordinary fermions, particularly for large 
£. In the other choice, one sets r s = 1. This presumably elevates the doublers 
well above the ordinary fermion modes, however, the optimal values of the 
parameters in this choice might receive O(moa s ) corrections, as suggested by 
tree-level and one-loop perturbative studies [8,11]. Therefore, if one takes the 
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choice of r s — 1, optimal values of the action parameters in principle must be 
tuned in a quark mass- dependent way. 

In this paper, we will discuss the tuning of the bare speed of light parame- 
ter v in a quenched calculation using tadpole-improved Wilson fermions on 
anisotropic lattices. The parameter v has to be tuned such that the pseudo- 
scalar meson energy-momentum dispersion relation reproduces its continuum 
form in the low- momentum limit. The dispersion relations of pseudo-scalar 
mesons are measured in our simulation for quark mass values ranging from the 
strange to the charm quark mass. The results of pseudo-scalar meson disper- 
sion relations at different values of v then enable us to extrapolate/interpolate 
to the optimized value of the bare speed of light parameter v for a given quark 
mass at a given gauge coupling (5. In order to measure the meson dispersion re- 
lations with better accuracy, asymmetric spatial lattice volumes are used which 
provide us with more non-degenerate (in the sense of energy) low-momentum 
modes. The quark action thus obtained can then be utilized in future studies 
on physical properties of hadrons with either light or heavy quarks. 

This paper is organized in the following manner. In Section 2, a particular form 
of clover-improved Wilson fermion action on anisotropic lattices is introduced. 
In Section 3, the calculation of the energy levels and dispersion relations for 
pseudo-scalar meson is discussed. This is performed for quark mass values 
ranging from the strange all the way to the charm quark mass at various values 
of gauge coupling and bare speed of light parameter v. By extrapolation or 
interpolation, the optimal values of the bare speed of light (denoted by u opt ) 
can then be determined for various values of (3 for a given bare quark mass 
parameter. In particular, we give the estimates for the optimal choice of v at 
the physical charm and strange quark mass values for a given (5. In Section 4, 
we will conclude with some general remarks. 



2 Improved Wilson Fermions on Anisotropic Lattices 

Consider a finite four-dimensional lattice with lattice spacing a M along the \x 
direction with \x = 0, 1, 2, 3. For definiteness, we denote a = a t and = a s for 
i — 1,2, 3. We will use £ = a s /a t to denote the bare aspect ratio of the asym- 
metric lattice. The quark actions on anisotropic lattices have been studied ex- 
tensively in the literature [8,9,10,27,14,11,28,23,16, 12,25,24,26]. Using these ac- 
tions, charmed meson spectrum [13,17], charmonium spectrum [14,16], charmed 
baryon spectrum [15,18] and hadron-hadron scattering lengths [19,20,21,22] 
have been studied. 

We start from the fermion action in the hopping parameter parametrization: 
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S = lf} X M Xy lf)y , 
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(2) 



Here we follow the notation as in Ref. [11], where we have made the choice 
r t = r s = 1 for the Wilson parameters. Another parameter ( = K s /K t is 
also commonly used in the literature. The forward and backward covariant 
derivatives on the lattice are given by: 



a fi V( l 'il) x = U fi (x)il> x+fi - ip x , 



(3) 



Using these definitions, one can rewrite the fermion action (2) in continuum- 
like notations: 



S = Y(ataM c) Mi c J4 c) , 
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where the the continuum fields and the bare quark mass m are given by: 



/2k 



m oa t = 1 - 3C ■ 



2k, ' 



For later convenience, we introduce the notation: 

1 P 

V = ^ ' 2k = 2/^ = m ° a * + ^ + 3z/ • 



(5) 



(6) 



We call the parameter v the bare speed of light parameter. The tuning of this 
parameter will be discussed in the remaining part of this paper using pseudo- 
scalar meson dispersion relations. Note that the critical bare quark parameter 
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depends explicitly on the parameter v even in the free case. This dependence 
also shows up qualitatively in our simulation. 



In quenched calculations, one usually needs to calculate the quark propagators 
at various valance quark masses. This amounts to different values of m or k 
for the same gauge field configuration. In this case, it is convenient to use the 
following fermion matrix: 



•^^xy ^xy^ + A X y 



A =S 

•^■xy u xy 



l/(2K max ) + p t ^(ToiFoi + Ps{0YlTvi + 0"23^23 + 0"31^3l) 

- H ^ [(! - nfp)Up(x)6 x+lt ,y + (1 + ln)Ul{x - n)6 x -^ y 



(7) 



where the coefficients are given by: 



* = 2 ' m= 2 ' a= 2k-2Z, 



Pt = v- 



(i + O 



v 

' ps= 2 



Here we have used the tree-level, zero quark mass relation: [11] 

cb = 1 , c E = —— , 



(8) 



(9) 



Note that, in principle the parameters cb and ce also have complicated depen- 
dence on the bare quark mass which we neglect in this study. In this notation, 
the bare quark mass dependence is singled out into parameter a and the ma- 
trix A remains unchanged when the bare quark mass is varied. Therefore, one 
could utilize the shifted structure of the matrix M. to solve for quark propaga- 
tors at various values of m (or equivalently k) at the cost of solving only one 
value of k = n rnaX) using the so-called Multi-Mass Minimal Residual (M 3 R for 
short) algorithm [29,30,31]. 

To implement the tadpole improvement, one replaces each spatial link Ui(x) 
by Ui{x)/u s while keeping the temporal link unchanged. 2 This results in the 
same fermion matrix (7) except that the parameters are replaced by: 



2 One can also introduce the tadpole improvement parameter ut for the temporal 
link. For large anisotropy £, this turns out to be irrelevant since the temporal lattice 
spacing is small enough and for all practical purposes, one can set ut = 1. 
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(10) 
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4u 2 s 
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2uj 



It is the quark action with these parameters that will be studied in this paper 
numerically. 



3 Simulation Results 

In this section, we present our numerical results for the study of the pseudo- 
scalar meson dispersion relations for various gauge coupling (3. Our main fo- 
cus lies upon the tuning of the bare speed of light parameter v for a given 
gauge coupling (3 and a given bare quark mass. The parameter v has to be 
tuned such that the lattice energy-momentum dispersion relations of pseudo- 
scalar mesons under investigation reproduce the continuum form in the low- 
momentum limit. To achieve this goal, one has to go through several proce- 
dures which we will describe in the following. 



3. 1 Simulation parameters and meson correlation functions 

The basic parameters of our simulation are summarized in Table 1. For the 
study of pseudo-scalar meson dispersion relations, it is advantageous to use 
lattices with asymmetric three- volume. This provides more non-degenerate (in 
the sense of its energy) low-momentum modes than the conventional symmet- 
ric volumes. All lattices in this study are of the size 6x9x12x50 except for 
the lattices at (3 = 3.0 where 8 • 12 • 16 ■ 50 lattices are studied. 3 To further 
check finite volume effects, a low statistics run (about 120 configurations) for 
(3 = 2.8 with larger lattice volumes was also performed. It turns out that the 
light meson mass values are somewhat modified but the final result of the 
optimal value of v remain compatible within errors (see Table 2). The aspect 
ratio is £ = a s /a t ~ £ = 5 for all lattices. The value of f3 ranges between 2.2 
and 3.0, roughly corresponding to spatial lattice spacing a s between 0.12 and 
0.27fm in physical units. For each particular value of (3, gauge field configura- 
tions are generated using the conventional pseudo-heatbath algorithms with 
over-relaxation. 



3 In our preliminary studies, 4 x 6 x 8 x 40 have also been simulated. We choose 
to present our results for larger lattices since they yield better accuracy for the 
pseudo-scalar meson dispersion relation measurements. 
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Table 1 



Simulation parameters for lattices (all with £ = 5) studied in this work. All lattices 
are of the size 6 x 9 x 12 x 50 except for j3 = 3.0 for which the lattice sizes are 
8 x 12 x 16 x 50. 





2.2 


2.4 


2.6 




r /a s = 1.76 


r /a s = 2.18 


r /a s = 2.48 


V 


0.86 


0.90 


0.94 


0.85 


0.90 


0.94 


0.87 


0.90 


0.93 


m cr (u) 


7.723(8) 


7.879(7) 


8.035(8) 


7.677(9) 


7.863(8) 


8.021(5) 


7.764(4) 


7.878(5) 


7.982(6) 


K 


0.0630 


0.0620 


0.0605 


0.0630 


0.0620 


0.0610 


0.0630 


0.0620 


0.0615 




0.0625 


0.0615 


0.0600 


0.0625 


0.0615 


0.0605 


0.0625 


0.0615 


0.0610 




0.0620 


0.0610 


0.0595 


0.0620 


0.0610 


0.0600 


0.0620 


0.0610 


0.0605 




0.0615 


0.0605 


0.0590 


0.0615 


0.0605 


0.0595 


0.0615 


0.0605 


0.0600 




0.0605 


0.0595 


0.0580 


0.0605 


0.0595 


0.0585 


0.0605 


0.0595 


0.0590 




0.0590 


0.0580 


0.0565 


0.0595 


0.0585 


0.0575 


0.0590 


0.0585 


0.0580 




0.0575 


0.0565 


0.0550 


0.0585 


0.0575 


0.0565 


0.0575 


0.0570 


0.0565 




0.0560 


0.0550 


0.0535 


0.0575 


0.0565 


0.0555 


0.0565 


0.0555 


0.0555 




0.0545 


0.0540 


0.0525 


0.0565 


0.0555 


0.0545 


0.0555 


0.0545 


0.0545 




0.0540 


0.0530 


0.0520 


0.0560 


0.0550 


0.0540 


0.0550 


0.0540 


0.0540 




0.0535 


0.0525 


0.0515 


0.0555 


0.0545 


0.0535 


0.0545 


0.0535 


0.0535 




0.0530 


0.0520 


0.0510 


0.0550 


0.0540 


0.0530 


0.0540 


0.0530 


0.0530 


(3 


2.8 


(6 ■ 9 ■ 12 ■ 


50) 


2.8( 


8- 12 ■ 16 ■ 


50) 


3.0 




r /a s = 3.13 


r /a s = 3.13 


r /a s = 4.13 


V 


0.87 


0.90 


0.93 


0.87 


0.90 


0.93 


0.87 


0.90 


0.93 


m cr (u) 


7.744(9) 


7.859(6) 


7.972(7) 


7.739(6) 


7.862(5) 


7.972(6) 


7.740(3) 


7.850(4) 


7.968(4) 


K 


0.0635 


0.0625 


0.0615 


0.0635 


0.0625 


0.0615 


0.0635 


0.0625 


0.0615 




0.0630 


0.0620 


0.0610 


0.0630 


0.0620 


0.0610 


0.0630 


0.0620 


0.0610 




0.0625 


0.0615 


0.0605 


0.0625 


0.0615 


0.0605 


0.0625 


0.0615 


0.0605 




0.0620 


0.0610 


0.0600 


0.0620 


0.0610 


0.0600 


0.0620 


0.0610 


0.0600 




0.0615 


0.0600 


0.0590 


0.0600 


0.0590 


0.0580 


0.0600 


0.0590 


0.0580 




0.0610 


0.0590 


0.0585 


0.0580 


0.0570 


0.0560 


0.0580 


0.0570 


0.0560 




0.0600 


0.0580 


0.0575 


0.0575 


0.0565 


0.0555 


0.0575 


0.0565 


0.0555 




0.0590 


0.0570 


0.0565 


0.0570 


0.0560 


0.0550 


0.0570 


0.0560 


0.0550 




0.0580 


0.0565 


0.0555 
















0.0570 


0.0560 


0.0550 
















0.0565 


0.0555 


0.0545 
















0.0560 


0.0550 


0.0540 















For gauge field configurations at a given value of (3, 3 different values of the 
bare speed of light parameter v as shown in Table 1 are studied. Quark prop- 
agators with zero and non-zero three-momenta are obtained using the M 3 R 
algorithm with wall sources for each dirac and color index. With the help of 
the M 3 R algorithm, by solving only one linear equation, one obtains the quark 
propagators with 12 (in the case of /3 — 3.0, only 8 values of k were taken) 
different values of k which correspond to different quark masses. The values of 
k are chosen such that the quark mass values range roughly from around the 
physical strange quark mass all the way up to the physical charm quark mass. 
The values of k for each parameter set ((3, v) are also tabulated in Table 1. 

In this paper, we focus on the single pseudo-scalar states with definite three- 
momentum. We define the pseudo-scalar and vector meson operators as fol- 
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lows: 

P(x,t) = g 1 (x,t) 75 g 2 (x,t), (11) 



where q±, q 2 (q~i, ?2)are quark field operators of two (possibly identical) flavors. 
Operators which create meson states with definite three-momentum k are then 
defined accordingly: 




where V 3 designates the three-volume of the lattice. Using these operators, 
one constructs the corresponding meson correlation function: 

C( PS \k,t) = (P(k,t)^P(k,0)) , (13) 



Using Wick's theorem, the above defined correlation function can be expressed 
in terms of the quark propagators: 



^ (P5) (k,t) = ^E^ 6 



V(2)ct6 
^aayt 



(14) 



where the Greek subscripts/superscripts in the solution vectors X and Y are 
Dirac indices while Roman subscripts/superscripts are color indices. The so- 
lution vectors X and Y are given by the inverse fermion matrix elements: 



ypb KA-l p-*x 

I f3ayt — Vt /3ayt;pbxO e ) 

x 

Xf3ayt = ^2 -^/3ayt;p6xO ■ (^) 
x 

The superscript (1) or (2) on these solution vectors indicates that the quark 
mass should be that of quark flavor 1 or 2 (possibly the same). These solution 
vectors are obtained by solving the linear equation of the fermion matrix M. 
with a suitable wall source. 



The energy values of a pseudo-scalar meson with definite three-momentum 
k (including zero-momentum) is obtained from their respective correlation 
functions C^ P5 ^(k, t) by finding the plateaus in their effective mass plots. In 
Fig. 1, we show the effective mass plots of a pseudo-scalar meson with zero 
three-momentum for /3 = 2.4, v = 0.85. The pseudo-scalar meson consists of a 
quark and an anti-quark with all possible bare quark mass value combinations 
(m 1 ,m 2 ). In Fig 1, we illustrate one of the situations with mi being fixed. 
The effective mass plots are shown for all m 2 values. Different lines in the plot 
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Fig. 1. The effective mass plots for the pseudo-scalar meson made up of a quark 
and an anti-quark for zero three-momenta at (3 = 2.4, v = 0.85. With one quark 
mass parameter m\ being fixed, different lines in the plot correspond to different 
values of the other quark mass m.2- The red horizontal bars in the plot indicates the 
ranges in which the mass of the meson are extracted. 



then correspond to different values of the other bare quark mass parameter 
m%. There are 12 lines in each of these windows which correspond to 12 dif- 
ferent values of the bare quark mass m 2 being simulated. It is seen that all 
effective mass plots develop nice plateaus at large temporal separation and ac- 
curate values of the pseudo-scalar energy E P s(mi,m 2 , k) ( and also the mass 
Mps(mi,m 2 ) = Eps(m,i, m 2 , k = 0) ) can thus be extracted. The red horizon- 
tal bars in the plot indicate the ranges in which the meson energy values are 
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extracted. The errors for the data points in this plot are analyzed using the 
standard jack-knife method. The intervals from which we extract the energy 
values are self-adjusted according to the minimum of \ 2 per degree of freedom. 
The quality of the effective mass plots for other parameter sets are similar. 



3.2 Obtaining the pseudo-scalar meson energy at fixed quark masses 



From the effective mass plots of pseudo-scalar meson correlation functions we 
obtain the energy values of a single pseudo-scalar meson with definite three- 
momentum k: Ep S (m 1: m 2: v, k). We thus have these energy values for each 
f3, u, k and all possible bare quark mass values (m 1 ,m 2 ) which we choose to 
calculate the meson correlation functions. Here the bare quark mass parameter 
is defined via: 

m =T. t. — -r^r • (16) 

2k 2k ct (v) v 1 



where K cr {v) is the critical hopping parameter at which the pion mass vanishes 
for a particular (/3,v). Note that this value depends on v for a given (3. For 
a given value of /3, the critical hopping parameter l/(2/t cr (z/) for each v is 
obtained by fitting the pion (made up of equal mass quarks) mass squared 
versus 1/{2k) using a quadratic function in the low quark mass region. From 
these fits, we obtain the critical value \/{2k ct {v)) for each v at a given (3. 

The reason that we choose the bare quark mass parameter m instead of the 
hopping parameter itself is the following. Our goal is to find the optimal values 
of v such that the pseudo-scalar meson exhibits the proper dispersion relation 
in the low-momentum region. Therefore, for a given value of /3, we want to fix 
the quark mass values and interpolate/extrapolate in v to obtain the optimal 
value of v at which the meson dispersion relation has the right form. This 
should be done for all possible quark mass values. It is better to perform this 
interpolation/extrapolation for fixed m instead of fixed hopping parameter 
pair k since the critical hopping parameters themselves depend explicitly on 
is, as is seen evidently from the tree-level relation Eq. (6). This dependence 
is also seen from our simulation. So for different values of u, the same value 
of k for different v in fact corresponds to different bare quark mass values. 
Therefore, it is more appropriate to interpolate/extrapolate in v for fixed bare 
quark mass parameter pair m. Note that the bare quark mass m as defined 
in (16) is an independent parameter of the quark action. In other words, no 
matter what the value of l/(2/t cr (z/)) comes out to be for each u, we could 
choose values of m, independent of u, since we are free to adjust the set of 
values for k. 
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Fig. 2. Interpolation of the pseudo-scalar meson energy squared E"p S {m\,m2, v i k) 
versus the quark mass parameter pairs {jn\^m%) to the common values (fh\,fh2) is 
shown for (3 = 2.4, v = 0.85. Only the zero three-momentum case is shown in this 
plot. Values of Ep S (mi,m2,i',k = 0) are shown as data points. Red lines are the 
quadratic interpolations around each value of fh using 6 points of m close to it. 



It turns out that our choices of the hopping parameters for different v are 
such that the range of m are roughly the same for a given (3 while the in- 
dividual values are not identical. Therefore, before we make any interpola- 
tion/extrapolation in v 1 which should be done at fixed m for all u, we first 
have to interpolate the energy values E PS (mi,m 2 ,v,]c) at different v to the 
the same quark mass values: (mi,m 2 ). In the analysis, we pick the values 
of fh to be the average of the three corresponding values at three different 
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v. Of course, the choice of the common values for fh at different v is some- 
what arbitrary and any other choice is equally well as long as the range of m 
roughly coincides with the ranges of the m for different v such that the inter- 
polation can be done reliably. The interpolation of the energy values squared 
Ep S (mi, mi-, v i k) is performed by a quadratic interpolation in the quark mass 
parameter using 4-6 neighboring points close to the values of fh. It is checked 
that all interpolations yield good fitting qualities. In Fig. 2, we show this in- 
terpolation for j3 = 2.4, v = 0.85. Values of Ep s (mi, m 2 , v, k) are shown as 
data points in the plot which correspond to zero three-momentum. The green 
points are the values of fh to which £p S (m 1 ,m 2 ,^, k) are interpolated. Red 
line segments are the quadratic interpolation in the quark mass parameter us- 
ing 4 — 6 points around each value of fh. The number of points being taken in 
each interpolation is determined by the condition of minimum \ 2 P er degree 
of freedom. The situation for non-zero momentum is similar. All interpolation 
yields good fitting quality. 

As the outcome of this procedure, we have all the quantities: Ep s (fhi,fh2, v, k) 
that are at the same sequence of points fh for different values of v. This proce- 
dure is performed for every value of (3 and for every three-momentum mode un- 
der investigation (see Eq. (refeq:ns) for the momentum modes being studied). 
In the discussion below, we will drop the bars on the quark mass parameters 
for simplicity with the understanding that all energy levels Eps(mi,m2, v, k) 
are already interpolated to the same set of (mi,m 2 ) for different v. 



3. 3 Extraction of Z parameter 

In this work, we utilize lattices with asymmetric volumes. This asymmetry 
helps to break the cubic symmetry in the momentum space and lifts the 
degeneracy of the meson energies. For example, by using a lattice of size 
6 x 9 x 12 x 50, we have 8 non-degenerate low-momentum modes, compared 
with only 4 with symmetric volumes. This technique also proves to be use- 
ful for the measurement of other momentum-dependent quantities like the 
hadron-hadron scattering phase shifts [32,33]. The three-momenta k in an 
asymmetric box of size L\ x L 2 x L 3 4 are quantized according to: 



with n = (711,712,713) G Z 3 being three-dimensional integers. In this work, the 
energy values of a meson with the following 8 three- momentum are measured: 



4 For definiteness, we pick L\ < L 2 < L$. 




(17) 
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n = (0, 0,0), (0,0,1), (0,1,0), (1,0,0), 
(0,1,1), (1,0,1), (1,1,0), (1,1,1). 



(18) 



For a given three-momentum k, the pseudo-scalar meson energy levels Ep s (mi,m 2 , 
are fitted according to the expected continuum dispersion relation: 

Ep s (m u m 2 ,v,\t) = M 2 PS (m l ,m 2 ,v) + Z PS (m 1 ,m 2 , v)k 2 , (19) 

in the low-momentum region. The fitting is performed for pseudo-scalar mesons 
with all possible bare quark mass combinations: (m 1; m 2 ). As a result of these 
fits, we obtain all the Z parameters of the pseudo-scalar mesons as a func- 
tion of two bare quark mass parameters: Zps(mi,m 2 , v). In particular, for the 
pseudo-scalar meson made up of quarks with the same mass, the correspond- 
ing Z parameter depends on one quark mass parameter, namely Zps(m, m, v). 
In Fig. 3, the linear fits of dispersion relations Ep s (m,m, u, k 2 ) versus k 2 for 
the pseudo-scalar meson with equal mass quark and anti-quark are shown for 
(3 = 2.4, v = 0.85. The straight lines in the plot represent the linear fits for 12 
values of bare quark parameter m. The linear fits utilize the low-momentum 
data points (including the zero-momentum point) according to Eq. (19) and 
the fitting range for each line is self- adjusted to yield the minimum x 2 for 
each degree of freedom. The slope of these lines then yield the parameters 
Zps(m, m, v) for all bare quark mass parameter m. Fitting qualities for other 
(/3, v) are quite similar. We have also tried another (conventional) way of ex- 
tracting the Z parameters, namely by using only the zero-momentum point 
and the lowest non- vanishing momentum point (n = (0,0,1) in this case). 
This is what has been done in the literature by many authors [12,24,25,26]. 
We find that the Z parameters are always better determined by using linear 
fits with more momentum points as compared with only two lowest momen- 
tum points. Therefore, we see the advantage of using asymmetric volumes for 
all our parameter sets. 

3.4 Finding the optimal values of v 

The optimal choice for the bare speed of light parameter v in the quark action 
is determined from the corresponding pseudo-scalar meson dispersion rela- 
tions, or more explicitly, from the Z parameters Zps(m,m, v) extracted from 
dispersion relations which is discussed the previous subsection. One requires 
that the dispersion relation of the pseudo-scalar meson made up of the same 
quark flavor reproduces its continuum counter-part in the low-momentum 
limit. That is to say, the optimal choice of v has to be such that: 



Zp S \m,m,v opt {m)\ = 1 . 



(20) 
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0.02 0.04 0.06 
k 2 

Fig. 3. Dispersion relations of the pseudo-scalar meson at (3 = 2.4, 0.85. The data 
points show the pseudo-scalar meson energy squared Ep s (m,rn,u,'k.) at a given 
three- momentum k versus k 2 . The straight lines are the linear fits to the data with 
the fitting range starting from the low-momentum end and self- adjusted according 
to the minimum x 2 P er degree of freedom. The slope of the lines then yield the 
desired Z parameters. Different lines in the plot indicate the linear fits for different 
bare quark mass parameter m. 

This yields the optimal value of v as a function of the bare quark mass pa- 
rameter: u op t(m). 

In practice, we use the values of Zps(m,m,i>) at different v and perform a 
linear extrapolation/interpolation in v for every value of m. This is shown in 
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Fig. 4. Determination of optimal speed of light parameter v ov t{yn) is shown for 
pseudo-scalar meson made of equal mass quark and anti-quark. Each small window 
corresponds to different values of bare quark mass parameter m. Data points in each 
window are the values of Zps(rn,m,u). They are fitted linearly versus v and the 
optimal values of v are determined by the condition: Zps[m, m, v pt(jn)] = 1 f° r each 
given m. The results of f op t(m) are shown as green points in each window together 
with the corresponding error bar. In each window (corresponding to different m), 
the quality of the linear fit is also indicated. 

Fig 4 in the case of j3 = 2.4. Each small window in this plot corresponds to 
different values of bare quark mass parameter m. Data points in each window 
are the values of Zps(m, m, v) obtained from the pseudo-scalar dispersion 
relations. They are fitted linearly versus v and the optimal values of v are 
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determined by the condition: Zps[m, m, u opt (m)] = 1. The results of v ovt {ra) 
for various to are shown as green points together with the corresponding error 
bars. In each window (corresponding to different to), the quality of the linear 
fit is also shown. 



3. 5 Finding the physical bare quark mass parameters 
Table 2 

Extracted physical bare quark mass parameters for the charm and the strange quark 
and the optimal values for the bare speed of light parameter v at physical quark 
masses for different (3. Two sets of data (I) and (II) for [3 = 2.8 corresponds to 
smaller (6 • 9 • 12 • 50) and larger (8 • 12 • 16 • 50) lattices, respectively. 






2.2 


2.4 


2.6 


2.8(7) 


2.8(77) 


3.0 


m ( t V) 


1.879(2) 


1.470(2) 


1.194(1) 


0.969(1) 


0.978(5) 


0.677(1) 


m ( F hy) 


0.217(3) 


0.163(1) 


0.155(1) 


0.120(2) 


0.15(2) 


0.073(2) 


v ov t(m? hV) ) 


0.68(3) 


0.71(3) 


0.754(5) 


0.75(1) 


0.74(3) 


0.84(1) 


Vo P t{m£ hy) ) 


0.940(4) 


0.933(3) 


0.928(3) 


0.98(2) 


0.91(4) 


0.99(1) 



For phenomenological reasons, one is particularly interested in the optimal 
parameters of v near bare quark mass values that correspond to the physical 
interesting cases. In particular, we are interested in the values of v opt (m) at 
physical bare strange quark and bare charm quark mass, namely at to = m^f^ 
or to = m!f hy \ To find out this correspondence, one has to investigate the 
bare quark mass dependence of the meson mass and fix the bare quark mass 
parameter to which corresponds to the physical case. 

To fix the physical bare quark mass parameters for the charm and the strange 
quark, we use the vector meson mass values. From the physical J/ ip mass and 
the physical $ meson mass, we can fix the bare quark mass parameters for 
the charm and the strange respectively. First, the vector meson mass squared 
at the optimal value of v: M v [m,m,i' opt (m)] is obtained by extrapolating to 
the optimal value of v. Once this quantity is at our disposal, we can perform 
extrapolation/interpolation in the bare quark mass to to find out the physical 
bare quark mass for the charm and the strange quark. We therefore extrapo- 
late My(m, m, v) versus Z PS (m, m, v) linearly for each given value of to. The 
situation is shown in Fig 5 for f3 = 2.4. Different small windows corresponds to 
different values of to. The linear extrapolations are shown by the straight lines 
in the windows. The extrapolated values of M v (m, m, v) at Zps(m, m,u) = 1, 
then give the quantities M v [to, to, u opt (m)]. The results of M v [m, m, u opt (m)] 
for all to are then utilized in the quark mass interpolation/extrapolation. 

In this work, we perform two quadratic fits for the meson mass versus the bare 
quark mass parameter to, one in the low quark mass region, the other in the 
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Fig. 5. Linear extrapolations of My(m,m,v) versus Zps(m,m,i/) for each value 
of m at (3 = 2.4. Different windows corresponds to different values of m. The linear 
extrapolations are shown by the straight lines in the windows. 

heavy quark mass region. We always take the fitting form to be: 

a^Mp S , v [m, m, u opt (m)] = A + Bm + Cm 2 . (21) 



In all our cases, we find the fit parameter A for the pseudo-scalar meson in the 
low quark mass region is always consistent with zero as it should be. In Fig. 6, 
we show this extrapolation for the vector meson mass at (3 = 2.4. The red line 
in the plot indicates a fit in the lower quark mass region. The green line is the 
corresponding fit in the heavy quark mass region. The corresponding fitting 
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Fig. 6. Quadratic extrapolations of My[m,m,u pt{i r n)] versus the bare quark mass 
m at (5 = 2.4. The red line indicates a fit in the lower quark mass region. The green 
line is the corresponding fit in the heavy quark mass region. The corresponding 
fitting ranges are also shown in the figure. The fitting ranges are self-adjusted to 
yield minimum x 2 P er degree of freedom. The blue horizontal line indicates the 
value for the physical J/tp meson while the pink horizontal line corresponds to the 
value of the physical <J? meson. The intersect with the green line and the red line 
then yields the estimate for m?' 9 ' and m^f hy \ respectively. 

ranges are also shown in the figure. The fitting ranges are self-adjusted to 
yield minimum \ 2 per degree of freedom. To obtain the physical charm quark 
mass parameter m!f hv " > , we draw a horizontal line in this figure at the physical 
J/ip mass: a 2 Mj,,. This is obtained by setting the scale using some physical 
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quantity. In this work, we choose the hadronic scale tq = 0.5fm (the so-called 
Sommer scale) to set the physical scale. 5 For different values of gauge coupling 
f3, the values of r /a s are known from the literature [2,3] which are also listed 
in Table 1 for reference. With this information, we know the physical meson 
masses in lattice unit. The blue horizontal line in the figure representing the 
value for physical J/ ip intersects with the green line and the intersection point 
then yields the estimate for the physical charm quark mass parameter m]? hy >. 
Similarly, the pink horizontal line is at the value of physical $ meson and the 
intersection point with the red line in the lower quark mass region yields the 
estimate for m < f hy \ The value of m!f hy) and m < f hy ^ thus obtained are listed in 
Table 2 for all (3. 

3. 6 Optimal values of v at physical quark mass parameters 

After fixing the physical bare quark mass parameters for both the charm 
and the strange, we can obtain the optimal values of the speed of light pa- 
rameter for these cases. In our notation, they correspond to v ovt (m!f hy) ) and 
^o- P t( m( f hy ^) , respectively. To get these values, we make interpolation/extrapolations 
of z / op*( m ) versus the bare quark mass parameter m using quadratic functions 
in these parameters. We choose the appropriate range (heavy quark mass re- 
gion and light quark mass region) for different cases. 

In Fig 7 we show the extrapolation of u opt {m) versus m in both the light and 
the heavy quark mass region at f3 — 2.4. The red line is the quadratic fit to the 
data points in the low quark mass region while the green line is the fit in the 
heavy quark mass region. The final result of u opt (m l f hy ^ ) and iy op t('m < f hv " > ) are 
indicated by the blue and the pink solid point with errors at the position of 
the physical charm and strange quark mass respectively. The results for other 
values of (5 are summarized in Table 2. 

For the (3 = 2.8 lattices, since the physical volume is somewhat small, one 
should check the size of finite size corrections. We therefore performed a low 
statistics run (about 120 gauge field configurations) for this (3. The final result 
is also tabulated in Table 2 labelled as f3 = 2.8(H). We see that the physical 
bare quark mass parameters changes, especially for the strange quark. How- 
ever, the optimal value for the parameter v at physical quark masses remain 
compatible within one or two standard deviations. Therefore, for the purpose 
of tuning the parameter z/, this result seems to indicate that the finite size 
corrections are not large. Of course, if one would use the action to calculate 
physical quantities and compare with the experimental values, it would be 
safer to use larger volumes, which is what we will do in the future. 

5 We have assumed that rg = 0.5fm exactly. Therefore, errors in this scale are not 
taken into account in the following error analysis. 
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Fig. 7. The result of Vopt^jn) is shown versus the bare quark mass m for (3 = 2.4. 
The red and the green line is a quadratic fit to the data points in the light and the 
heavy quark mass region respectively. At the position of the physical bare charm 
quark mass mf hv ' and the physical bare strange quark mass m^f hv \ the result of 
v pt{ m c' hv " > ) an d Vopt{ms° hy ^) are shown by the blue and the pink solid point with 
errors, respectively. 



4 Conclusions 



In this paper, we present a systematic numerical analysis on the tuning of the 
bare speed of light parameter v in the tadpole improved anisotropic Wilson 
quark action. The tuning is done in a quark mass dependent way with quark 
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mass values ranging from the strange to the charm. The optimal values of v 
are obtained for various values of (3 using the pseudo-scalar meson dispersion 
relations. With the help of the anisotropic lattices with asymmetric volumes, 
the dispersion relations can be measured with good accuracy. Using the tad- 
pole improved anisotropic Wilson action with these optimized parameters, a 
quenched calculation can then be performed to study properties of hadrons 
made up of either light or heavy quarks. Therefore, with the same improved 
quark action one can study hadron spectrum and other physical properties 
in a wide range of quark masses. We hope to come to this issue in the near 
future. 

Acknowledgments 

We would like to thank Prof. H. Q. Zheng and Prof. S. L. Zhu of Peking 
University for helpful discussions. Our thanks also goes to Dr. J. P. Ma at 
Institute of Theoretical Physics, Academia Sinica and Dr. Y. Chen at Institute 
of High Energy Physics, Academia Sinica for their stimulating discussions. 

References 

[I] C. Morningstar and M. Peardon. Phys. Rev. D, 56:4043, 1997. 
[2] C. Morningstar and M. Peardon. Phys. Rev. D, 60:034509, 1999. 
[3] C. Liu. Chinese Physics Letter, 18:187, 2001. 

[4] C. Liu. Communications in Theoretical Physics, 35:288, 2001. 

[5] C. Liu. In Proceedings of International Workshop on Nonperturbative Methods 
and Lattice QCD, page 57. World Scientific, Singapore, 2001. 

[6] C. Liu and J. P. Ma. In Proceedings of International Workshop on 
Nonperturbative Methods and Lattice QCD, page 65. World Scientific, 
Singapore, 2001. 

[7] C. Liu. Nucl. Phys. (Proc. Suppl.) B, 94:255, 2001. 

[8] A.X. El-Khadra, A.S. Kronfeld, and P.B. Mackenzie. Phys. Rev. D, 55:3933, 
1997. 

[9] M. Alford, T. R. Klassen, and P. Lepage. Phys. Rev. D, 58:034503, 1998. 
[10] T. R. Klassen. Nucl. Phys. (Proc. Suppl.) B, 73:918, 1999. 

[II] J. Harada, A.S. Kronfeld, H. Matsufuru, N. Nakajima, and T. Onogi. Phys. 
Rev. D, 64:074501, 2001. 



21 



[12] J. Harada, H. Matsufuru, T. Onogi, and A. Sugita. Phys. Rev. D, 66:014509, 
2002. 

[13] P.B. Mackenzie, S.M. Ryan, and J.N. Simone. Nucl. Phys. (Proc. Suppl.) B, 
63:305, 1998. 

[14] P. Chen. Phys. Rev. D, 64:034509, 2001. 

[15] R. Lewis, N. Mathur, and R.M. Woloshyn. Phys. Rev. D, 64:094509, 2001. 
[16] M. Okamoto et al. Phys. Rev. D, 65:094508, 2002. 
[17] A. Juettner and J. Rolf. Phys. Lett. B, 560:59, 2003. 

[18] Y. Nemoto, N. Nakajima, H. Matsufuru, and H. Suganuma. Phys. Rev. D, 
68:094505, 2003. 

[19] C. Liu, J. Zhang, Y. Chen, and J.P. Ma. Nucl. Phys. B, 624:360, 2002. 

[20] G. Meng, C. Miao, X. Du, and C. Liu. hep-lat/0309048, 2003. 

[21] C. Miao, X. Du, G. Meng, and C. Liu. hep-lat/0403028, 2004. 

[22] X. Du, C. Miao, G. Meng, and C. Liu. hep-lat/0404017, 2004. 

[23] Junhua Zhang and C. Liu. Mod. Phys. Lett. A, 16:1841, 2001. 

[24] T. Umeda et al. Phys. Rev. D, 68:034503, 2003. 

[25] S. Hashimoto and M. Okamoto. Phys. Rev. D, 67:114503, 2003. 

[26] Justin Foley, Alan O Cais, Mike Peardon, and Sinead M. Ryan, hep-lat/0405030, 
2004. 

[27] Junko Shigemitsu Stefan Groote. Phys. Rev. D, 62:014508, 2000. 

[28] H. Matsufuru, T. Onogi, and T. Umeda. Phys. Rev. D, 64:114503, 2001. 

[29] A. Frommer, S. Giisken, T. Lippert, B. Nockel, and K. Schilling. Int. J. Mod. 
Phys. C, 6:627, 1995. 

[30] U. Glaessner, S. Guesken, T. Lippert, G. Ritzenhoefer, K. Schilling, and 
A. Frommer. hep-lat/9605008. 

[31] B. Jegerlehner. hep-lat/9612014. 

[32] X. Li and C. Liu. Phys. Lett. B, 587:100, 2004. 

[33] X. Feng, X. Li, and C. Liu. Phys. Rev. D, 70:014505, 2004. 



22 



